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Abstract 



We determine the CKM matrix element \V c b\ using a sample of 3.33 million 
BB events in the CLEO detector at CESR. We determine the yield of recon- 
structed B — > D* + lv decays as a function of w = vb • «d*, and from this we 
obtain the differential decay rate dF/dw. By extrapolating dT/dw to w = 1, 
the kinematic point at which the D* is at rest relative to the B, we extract the 
product |V^b|.F(l), where F(l) is the form factor at w = 1 and is predicted ac- 
curately by theory. We find F(l)\V cb \ = 0.0424± 0.0018(stat.) ±0.0019(syst.). 
We also integrate the differential decay rate over w to obtain B(B — ► D* + iu) = 
(5.66 db 0.29 db 0.33)%. All results are preliminary. 
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I. INTRODUCTION 



The CKM matrix element \V c b\ sets the length of the base of the famous unitarity trian- 
gle. One strategy for determining |V^,| uses the decay B — > D*iv. The rate for this decay, 
however, depends not only on \V c b\ and well-known weak decay physics, but also on strong 
interaction effects, which are parametrized by form factors. These effects are notoriously 
difficult to quantify, but Heavy Quark Effective Theory (HQET) offers a method for calcu- 
lating them at the kinematic point at which the final state D* is at rest with respect to the 
initial B meson (w = vd* ■ v B = 1, and is the relativistic boost 7 of the D* in the B rest 
frame). In this analysis, we take advantage of this information: we measure dT/dw for these 
decays, and extrapolate to obtain the rate at w = 1. The rate at this point is proportional 
to [|V^b|F(l)] 2 where F(w) is the form factor. Combined with the theoretical results, this 
gives \V cb \. 

The analysis uses B° — > D* + iv decays and their charge conjugates (charge conjugates 
are implied throughout this paper). We divide the reconstructed events into bins of w. In 
each bin we extract the yield of D*£u decays using a fit to the distribution cos6b-d*£, where 



cos6 B _ DH = —2- (1) 

2|Pb||Pd*£| 



The angle cos8b-d*£ is thus the reconstructed angle between the D*-lepton combination 
and the B meson, computed with the assumption that the only missing mass is that of the 
neutrino. This distribution distinguishes B — > D*lv decays from decays such as B — > D**lv, 
since D*iv decays are concentrated in the physical region, —1 < cos6b-d*£ < 1, while the 
larger missing mass of the D**iv decays allows them to populate cos8b-d*£ < — 1- Given 
the D*lv yields as a function of w, we fit for a parameter describing the form factor and the 
normalization at w — 1. This normalization is proportional to the product [| V^|-F(l)] 2 . 



II. EVENT SAMPLES 

We do our analysis with 3.33 million BB events (3.1 fb _1 ) produced on the T(45) 
resonance at the Cornell Electron Storage Ring and detected in the CLEO II detector ||T| . In 
addition, the analysis uses a sample of 1.6 fb _1 of data collected slightly below the T(4S") 
resonance for the purpose of subtracting continuum backgrounds. 

The analysis uses events from a GEANT-based [|| Monte Carlo simulation to provide the 
distribution of D*lv events in cos8b-d*£ an d to provide information on some backgrounds. 
In the simulation, D*lv decays are modeled using a linear form factor (for hA 1 {w)) with the 
parameters measured in a previous CLEO analysis ||. The signal includes events with final 
state radiation as modeled by PHOTOS jlj. We simulate other form factors by reweighting 
this sample. Non-resonant B — ► D*irlv decays are modeled using the results of Goity and 
Roberts [Q, and B — > D**lv decays are modeled using the ISGW2 form factors. In 
the following, we refer to B — ► D**£u and non-resonant B — > D*tt£u collectively as D*Xiv 
decays. 
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III. EVENT RECONSTRUCTION 



We use the decay chain D* + — ► D G -n + followed by D° — ► fT~7r + . We first combine kaon 
and pion candidates in hadronic events to form D° candidates. Signal events lie in the 
window \itik-k ~ 1 .865 1 < 0.020 GeV. We then add a slow ir + to the D° candidate to form 
a D* + . The K and tc are fit to a common vertex and then the slow n and D are fit to a 
second vertex using the beam spot constraint. This second vertex constraint improves the 
resolution in Am = ttik-k^ ~ m Kn by about 20%. We require \Am — 0.14544| < 0.002 GeV. 

Electrons are identified using the ratio of their energy deposition in the Csl calorimeter 
to the reconstructed track momentum, the shape of the shower in the calorimeter, and 
their specific ionization in the drift chamber. Our candidates lie in the momentum range 
0.8 < p e < 2.4 GeV. Muon candidates penetrate two layers of steel in the solenoid return 
yoke (~ 5 interaction lengths). Only muons above about 1.4 GeV satisfy this requirement; 
we therefore demand that they lie in the momentum range 1.4 < p^ < 2.4 GeV. The charge 
of the lepton must match the charge of the kaon, and be opposite that of the slow pion. 

Exact reconstruction of w requires knowledge of the flight direction of the B meson. 
While this is unknown, our knowledge of cos 9b~d*£ limits it relative to that of the D* — t 
combination. We therefore compute w using the directions at each end of the range and we 
then take the average. The typical resolution in w is 0.03. We divide our sample into 10 
equal bins from 1.0 to 1.51, where the upper bound is just beyond the kinematic limit of 
1.504. In a few events, the reconstructed w falls outside its kinematic range; we assign these 
to the first or last bin as appropriate. In the high w bins, we suppress background with no 
loss of signal efficiency by restricting the angle between the D* and the lepton. 

IV. EXTRACTING THE D*lv YIELDS 
A. Method 

At this stage, our sample of candidates contains not only D*£u events, but also B — > 
D*Xiv decays and various backgrounds. In order to disentangle the D*iv from the D*Xiv 
decays, we use a binned maximum likelihood fit |JJ to the cosQb-d*i distribution. In this 
fit, the normalizations of the various background distributions are fixed and we allow the 
normalizations of the D*lv and the D*Xiv events to float, with the constraint that both 
normalization be positive (or zero) and that the total event yield matches that of the data. 

The distributions of the D*lv and D*X£u decays come from our signal Monte Carlo. 
The backgrounds, and how we obtain their cos 6 b-d*( distributions and normalizations, are 
described in the next section. 

B. Backgrounds 

There are several sources of events other than B — > D*lv and B — > D*X£u. We divide 
these backgrounds into five classes: continuum, combinatoric, uncorrelated, correlated and 
fake lepton. 
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FIG. 1. The Am distribution of events in the fifth w bin, for all cos6b-d*£- The data (solid 
squares) are superimposed with the combinatoric background distribution (dashed curve) and the 
sum of the background and the D* + signal (solid histogram). The cos9b-d*£ distribution of the 
combinatoric events is taken from events in the region 0.155 < Am < 0.165 GeV. 

1. Continuum Background 

At the T(4S f ) we detect not only resonance events, but also non- resonant events such 
as e + e~ — > qq. This background contributes about 4% of the events within the range 
— 1 < cos9b-d*£ < 1 (the "signal region"). In order to subtract background from this source, 
CESR runs one-third of the time slightly below the T(4S) resonance. For this continuum 
background, we use the cos9b-d*£ distribution of events in the off-resonance data scaled by 
the ratio of luminosities and corrected for the small difference in the cross sections at the 
two center of mass energies. 

2. Combinatoric Background 

Combinatoric background events are those in which one or more of the particles in the 
D* candidate does not come from a true D* decay. This background contributes 6% of the 
events in the signal region. 

We take the cos8b-d*£ distribution of combinatoric background events from the high 
Am sideband (0.155 < Am < 0.165 GeV). Their normalization comes from a fit to the Am 
distribution in which we assume a background shape of the form n(Am — m 7r ) 1_c , and vary 
n, c, and the normalization of the signal peak. The lineshape for the peak is taken from 
Monte Carlo. This fit is shown for a representative w bin in Figure [I]. 
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3. Uncorrelated background 



Uncorrelated background, which accounts for approximately 4% of the events in the 
signal region, arises when the D* and lepton come from the decays of different B mesons. 
Most of this background consists of a D* meson combined with a secondary lepton (that is, 
a lepton from the chain b — > c — > slv) because primary leptons from the other B have the 
wrong charge correlation. Uncorrelated background events can also arise, however, when the 
B° and B° mix or when a D* from the upper- vertex (that is, from the c in the decay chain 
b — > ccs) is combined with a primary lepton. 

We obtain the cos 9b-d*v distribution of this background by simulating each of the vari- 
ous sources and normalizing each one appropriately. To normalize, we use the inclusive D* 
production rate observed in our data, the measured primary and secondary lepton decay 
rates |J, the estimated decay rate for modes in the B — > D^D*K^ family, and the mea- 
sured B° — B° mixing rate [ 10 1 . Since the cos6b-d*£ distribution depends somewhat on the 



momentum of the D*, we normalize the D* sources separately in low and high momentum 
bins. 



4- Correlated Background 

Correlated background events are those in which the D* and lepton are daughters of the 
same B, but the decay was not B — > D*£u or B — ► D*Xlv. The most common sources are 
B — > D*tv followed by leptonic r decay, and B — ► D*D S followed by semileptonic decay of 
the D s . This background accounts for less than 0.5% of the events in the signal region and 
is provided by Monte Carlo simulation. 



5. Fake Lepton Background 

Fake lepton background arises when a hadron is misidentified as a lepton and is then 
used in our reconstruction. A preliminary study indicates that this background is small, and 
we ignore it. 



C. B -» D*lv Yields 

Having obtained the distributions in cos Ob-dh of the signal and background components, 
we fit for the yield of D* + iu events in each w bin. Two representative fits are shown in 
Figure 0. The quality of the fits is good, as is agreement between the data and fit distributions 
outside the fitting region. We summarize the D*£u and D*Xiv yields in Figure |^. 

In order to test the quality of our fit and modeling of the signal and backgrounds, we 
plot a variety of distributions. Figure f| shows the D* energy distribution and the lepton 
momentum spectrum. We find good agreement between the data and our expectations. 
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FIG. 2. The event yields (circles) in the third and ninth uu bins with the results of the fit 
superimposed. The fit range is —8 < cos6b-d*£ < 1-5, and is indicated with the solid circles. 
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FIG. 3. The D*tv and D* + X£u yields in each w bin. 
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FIG. 4. The energy distribution of D* + candidates (left) and the lepton momentum spectrum 
for D* + H.v candidates (right) in the signal region for all w bins combined. 



V. THE \V cb \ FIT 
The partial width for B — > decays is given by 



dT G\ 



dw 487r 3 



(m B - m D *) 2 m 3 D *Vw 2 -l(w + l) 2 ( 1 + 4 



w 



w + 1 



1 — 2wr + r 
(1-r) 2 



(2) 



where and m^* are the S and D* meson masses, r = ttid* /uiq, and the form factor 
F(w) is given by 



F(w) 



Hl + Hj + H 2 

^ 1+4^J i (1— r) 2 



(3) 



The ifj are the helicity form factors and are given by 

w — 1 



Hjw) 



1 — r 



Vl — 2ur + r 2 / iw — 1 

1 - W r-Ri(w) 



1 — r 

Vl — 2u>r + r 2 
1 — r 



w + 1 



/iw - 1 
w + 1 



RiH 



(4) 
(5) 

(6) 



The form factor h Al (w) and the form factor ratios Ri(w) = hy{w)/h J s L1 {w) and R2(w) = 
{Iia a {w) + rhA 2 (w)) /h Al {w) have been studied both experimentally and theoretically. A 
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CLEO analysis H measured these form factors under the assumptions that h^iw) is linear 
as a function of w and that R\ and R 2 are independent of w. CLEO found 

i?i = 1.18 ±0.30 ±0.12, (7) 
R 2 = 0.71 ± 0.22 ± 0.07 and (8) 

- d ^(w = 1) =p 2 = 0.91 ±0.15 ±0.06 (9) 

dw 

with the correlation coefficients C{p 2 , R x ) = 0.60, C{p 2 , R 2 ) = -0.80 and C(R U R 2 ) = -0.82. 

i?i(l) and R 2 {1) have been computed using HQET and QCD sum rules with the results 
i?i(l) = 1.27 and #2(1) = 0.8 and estimated errors of 0.1 and 0.2 respectively [jTTf] , in good 
agreement with the experimental results. Ri(w) and R 2 (w) are expected to vary weakly 
with w. Most importantly for this analysis, F(l)(= /^(l)) is relatively well-known, thereby 
allowing us to disentangle it from |V cfc |. We will use F(l) = 0.913 ± 0.042 Q. 

Recently, dispersion relations have been used to constrain the shapes of the form factors 
TB| , [|14| . Rather than expand the form factor in w, these analyses expand in the variable 



z = ( V /^7TI_ y/2)/(y/wTT+ V2). Ref. [[TJ obtains the results: 

h Al {w) = 1 - 8p 2 z + (53p 2 - I5)z 2 - (231p 2 - 91)z 3 (10) 

i?iH = i?i(i) - 0.120 - 1) + 0.05O - 1) 2 (11) 

R 2 (w) = R 2 (l) + 0.11(w - 1) - 0.06(w - l) 2 . (12) 

In our analysis, we assume that the form factor has the functional form derived from 
dispersion relations given in Equations |10|, |Tl| and [l|. We fit our yields as a function of w for 
F(l)|Vtf,| and p 2 , keeping -Ri(l) and ^2(1) fixed at their measured values. Our fit minimizes 

10 \T\Tobs _ y-10 f ..AT.]2 

X 2 = E- 2 , (13) 

i=l ^j^obs 

where N° bs is the yield in the i th w bin, Nj is the number of decays in the j th w bin, and 
the matrix ey accounts for the reconstruction efficiency and the smearing in w. Explicitly, 

r dT 

Nj = Afo N Ti , S )B(D* + -> D°7T + )B(D° -> K-ti + )t b0 / dw— (14) 

Jwj dw 

where t b o is the B° lifetime JTJ], B(D* + — > D°7r+) is the D* + — > D°7r + branching frac- 
tion i3(-D — > Kit) is the D° — > i^~7r + branching fraction iV-f^s) is the number of 
T(AS) events in the sample, and /oo represents the T(45) — > B°B° branching fraction [16|. 
We use the result of Ref. |16| for / o as a constraint in the fit. 
The result of the fit is shown in Figure |5|. We find 

\V cb \F(l) = 0.0424 ± 0.0018 (15) 
p 2 = 1.67 ±0.11 and (16) 
X 2 = 3.1/8 dof. (17) 

with a correlation coefficient between |Vq&|.F(1) and p 2 of 0.90. These parameters give T = 
0.0366 ± 0.0018 ps _1 . The quality of the fit is excellent. We note that the slope is higher 
than that found in the previous CLEO analysis |17| because of the curvature introduced into 
our form factor. If we use a linear form factor and the same subset of the data, we obtain 
results compatible with the earlier analysis. 
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FIG. 5. The results of the fit to the w distribution. The upper figure shows the yields (solid 
circles) with the results of the fit superimposed (histogram). The lower figure displays \V c b\F(w), 
where the data points (solid circles) are derived from the yields after correcting for efficiency, 
smearing, and all terms in the differential decay rate apart from \V c b\F(w). The curve shows the 
result of the fit. 
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VI. SYSTEMATIC UNCERTAINTIES 



The systematic uncertainties are summarized in Table |. The dominant systematic un- 
certainties arise from our background estimations and from our knowledge of the slow pion 
reconstruction efficiency 

A. Background uncertainties 

We test our procedure for subtracting combinatoric background by applying it to Monte 
Carlo simulated events. We assign a systematic error based on the difference between the 
results obtained using the "true" background and those obtained using the same procedure 
that we apply to our data. We also include the statistical uncertainty of the study. 

The main source of uncertainty from the uncorrelated background is the normalization 
of the various contributions. Of these, the most important is the branching fraction of the 
B -> D^D*K^ decays, which we vary by 50%. Smaller effects arise from the primary and 
secondary lepton rates and from the uncertainty in B° — B° mixing. 

We assess the uncertainty arising from the correlated background by varying the branch- 
ing fractions of the contributing modes. 

B. Slow 7r reconstruction uncertainty 

A major source of uncertainty for the analysis is the efficiency for reconstructing the slow 
pion from the D* decay. This efficiency is low near w = 1 and increases rapidly over the 
next few w bins. We have explored the efficiency as a function of the event environment 
and as a function of hit resolution, hit efficiency, material, and the charge division resolution 
in one of the inner drift chambers. The last of these is important because very soft tracks 
often make use of the charge division information for reconstruction of the z component 
(beam direction) of their trajectory. The uncertainty in |V^,|.F(1) and the decay width are 
dominated by uncertainties in the amount of material in the inner detector (2.3%) and the 
drift chamber hit efficiency (0.8%). 

C. Other uncertainties 

The efficiency for identifying electrons has been evaluated using radiative bhabha events 
embedded in hadronic events, and has an uncertainty of 2.4%. Similarly, the muon identifi- 
cation efficiency has been evaluated using radiative mu-pair events, and has an uncertainty 
of 1.4%. The total uncertainty in lepton identification, weighted by the electron and muon 
populations, is 2.1%. Separate electron and muon analyses of our data give consistent results. 

The B° momentum is measured directly in the data using fully reconstructed hadronic 
decays, and is known on average with a precision of 0.0016 MeV. Variation of the momentum 
in our reconstruction slightly alters the cos8b-d*£ distribution that we expect for our signal, 
and it therefore changes the yields obtained from the cos 9b-d*£ fits. Likewise, CLEO has 
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measured the B° meson mass |T8[ and when we vary the mass within its measurement error, 
we find a small effect on the yields. 

We determine the tracking efficiency uncertainties for the lepton and the K and 7r forming 
the D° in the same study used for the slow pion from the D* + decay. These uncertainties 
are confirmed in a study of 1-prong versus 3-prong r decays. 

Finally, our analysis requires that we know the cos 9 b-d*i distribution of the D*X£u con- 
tribution. This distribution in turn depends on both the branching fractions of contributing 
modes and on their form factors. Variation of all of these branching fractions and form 
factors is not only cumbersome, but out of reach given the poor current knowledge of these 
modes. Instead, we note that the B — ► D*tt£u and B —>■ D\lv modes are the ones with the 
most extreme cos 8 'b-d* t distributions (the largest mean and the smallest). These distribu- 
tions are shown in Figure |6|. We therefore repeat the analysis, first using pure B — > D*ixtv 
to describe our D*X£v decays and then using pure B — > D x tv to describe these decays, and 
we take the larger of the two excursions as our systematic error. 



D. Sensitivity to R x {\) and R 2 (l) 

The form factor ratios -Ri(l) and #2(1) affect the lepton spectrum and therefore the frac- 
tion of events satisfying our 0.8 GeV electron and 1.4 GeV muon momentum requirements. 
To assess this effect, we vary -R(l) and -R(2) within their measurement errors, taking into 
account the correlation between them. 
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TABLE I. The fractional systematic uncertainties. 
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VII. CONCLUSIONS 

We have fit the w distribution of B — > D*£z/ decays for the slope of the form factor and 
\V cb \F(l). We find 

\V cb \F(l) = 0.0424 ± 0.0018 ± 0.0019 (18) 
p 2 = 1.67 ±0.11 ±0.22 (19) 

(20) 

with a correlation coefficient between |Vd,|.F(l) and p 2 of 0.90. These parameters imply the 
decay rate 

r = 0.0366 ±0.0018 ±0.0023 ps _1 . (21) 
and the branching fraction 

B(B -> D* + £v) = (5.66 ± 0.29 ± 0.33)%. (22) 

Our result implies 

\V cb \ = 0.0464 ± 0.0020(stat) ± 0.0021 (syst.) ± 0.0021 (thear.), (23) 

where we have used F(l) = 0.913±0.042 fl2|| . This is consistent with previous measurements 
of \V cb \, but is somewhat higher. The analysis benefits from small backgrounds and good 
resolution in w. These results are preliminary. 
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